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X-ray spectroscopy is a sensitive probe of stellar winds. X-rays originate from optically thin 
shock-heated plasma deep inside the wind and propagate outwards throughout absorbing cool 
material. Recent analyses of the line ratios from He-like ions in the X-ray spectra of 0-stars 
highlighted problems with this general paradigm: the measured line ratios of highest ions are 
consistent with the location of the hottest X-ray emitting plasma very close to the base of 
the wind, perhaps indicating the presence of a corona, while measurements from lower ions 
conform with the wind-embedded shock model. Generally, to correctly model the emerging X- 
ray spectra, a detailed knowledge of the cool wind opacities based on stellar atmosphere models 
is prerequisite. A nearly grey stellar wind opacity for the X-rays is deduced from the analyses 
of high-resolution X-ray spectra. This indicates that the stellar winds are strongly clumped. 
Furthermore, the nearly symmetric shape of X-ray emission line profiles can be explained if 
the wind clumps are radially compressed. In massive binaries the orbital variations of X-ray 
emission allow to probe the opacity of the stellar wind; results support the picture of strong 
wind clumping. In high-mass X-ray binaries, the stochastic X-ray variability and the extend of 
the stellar-wind part photoionized by X-rays provide further strong evidence that stellar winds 
consist of dense clumps. 



1 Introduction 

All OB-stars with spectral types earlier than B1.5 
have X-ray luminosities Lx exceeding lO'^"'^ ergs~^. 
Remarkably, their X-ray luminosity roughly cor- 
relates with stellar bolometric luminosities as 
Lx~10~''iboi (Seward et al. 1979, Berghofer et al. 
1994). The high-resolution X-ray spectra of two O- 
stars, CPup and COri, are shown in Figure 1. The 
spectra are dominated by strong emission lines that 
reflect stellar wind abundances, e.g. nitrogen lines 
are more prominent in the spectrum of the more 
evolved star ^ Pup. X-ray spectra of single non- 
magnetic massive stars can be explained as emit- 
ted by optically thin plasma heated to temperatures 
fcTx<0.6keV (e.g. Wojdowski & Schulz 2005). 

There is a number of hypotheses to explain how 
X-rays are generated in stellar winds. Lucy & White 
(1980) put forward a theory where X-rays originate 
from bow shocks around a population of blobs (see 
also Cassinelli et al., these proceedings). Cassinelli 
& Swank (1983) proposed an analogy with coronal 
stars where X-rays originate from magnetically con- 
fined regions at the wind base. Owocki et al. (1988) 
performed numerical simulations of the evolution of 
instabilities in radiatively driven stellar winds and 
predicted a development of reverse shocks with post- 
shock temperatures in the range 10^ — 10'' K (i.e. 
as inferred from the observed X-ray spectra). Feld- 
meier et al. (1997) developed a model explaining 
the X-rays as originating from cooling zones behind 
shock fronts and were able to reproduce the low- 
resolution Rosat spectra of Q Pup and C Ori. 



In massive binaries a copious amount of X-ray 
emission is produced by the collision of two pow- 
erful stellar winds. X-ray spectra of colliding wind 
binaries reveal that the X-ray emitting plasma de- 
parts from coUisional equilibrium (e.g. Pollock et al. 
2005). Colliding wind systems tend to have some- 
what higher X-ray luminosities and their X-ray spec- 
tra can be "harder" compared to the single stars. 
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and coUisional terms are omitted). Neglecting the 
limb darkening, the average intensity can be sub- 
stituted by the diluted photospheric flux, W{r)H^. 
If the photospheric flux H^, is known, the dilution 
factor W can be inferred from TZ. This constrains 
the radius i?fii. where the X-ray emitting plasma is 
located. 



Figure 1: High- resolution X-ray spectra of CPup 
(top panel) and ^Ori (bottom panel) ob- 
tained with XMM-Newton. Prominent 
emission lines are identified. 

Another type of X-ray emitting plasma is found 
in high-mass X-ray binary systems (HMXBs) that 
consist of an OB supergiant and a degenerate com- 
panion, usually a neutron star (NS) deeply embed- 
ded in the stellar wind of the primary. Stellar wind 
accretion onto a NS powers strong X-ray emission 
that has a characteristic power-law spectrum, while 
the surrounding stellar wind is photoionizcd by this 
X-ray emission. 

Independent on the formation mechanism, the X- 
rays registered by the observer are attenuated in the 
stellar wind and interstellar medium. The trans- 
port of X-rays in stellar winds can often be treated 
as pure absorption (Baum et al. 1992, Hillier et al. 
1993). This leads to a great simplification of the ra- 
diative transfer problem and allows modeling even 
in the case when the absorbing medium is inhomo- 
geneous. 



For wavelengths in the observable part of the UV 
the flux can be, in principle, directly inferred from 
observations. Unfortunately, especially for giant and 
supergiant stars, the photospheric fluxes are often 
contaminated by wind lines. E.g. Ovi resonance 
doublet coincides with the Af^i transition at 1033 A 
for the Mgxi line. The Ovi doublet itself can only 
be reproduced by stellar atmosphere models if an X- 
ray field causing Auger ionization is included (e.g. 
Oskinova et al. 2006). This example highlights the 
need of a full radiative transfer treatment for the 
correct interpretation of the fir line ratios - a task 
that has not been accomplished yet. 



2 X-ray emission of single stars as 
a probe of clumped stellar wind 

2.1 Line ratios in spectra of He-like ions 

One of the key spectral diagnostics is provided by the 
ratio of line fluxes from He- like ions. These ions show 
characteristic '/ir triplets' of a forbidden (f), an inter- 
combination (i) and a resonance (r) line. In OB-type 
stars, the TZ =f /i line ratio is sensitive to the local 
mean intensity of the UV radiative field, J^. The UV 
photons with wavelength Af_,i excite the metastable 
^S level to the level, so that TZ(x{l + (/)(J^))-\ 
where (p^Ju) is the photoexcitation rate (constants 



Commonly, the fluxes at the wavelengths of in- 
terest are obtained from stellar atmosphere models. 
Leutenegger et al. (2006) analyzed the lines of He- 
like ions in the spectra of four 0-type giants and su- 
pergiants. The parameters of the 0-stars were taken 
from Lamers & Leitherer (1993) and used to select a 
stellar atmosphere model from the TLUSTY static 
plane-parallel model grid. Leutenegger et al. find 
that the minimum radius of X-ray formation is typ- 
ically in the range of 1.25< i?fir/-R*< 1.67. 
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Figure 2: Logarithm of the radius were the ra- 
dial optical depth of the wind becomes 
unity in CPup according to PoWR stel- 
lar atmosphere models. Thick line: 
model with M=2.5x 10"^ Afo/yr; thin 
line: M=8.7x 10"^ Mg/yr. The radii of 
the X-ray emission derived from different 
He-like ions arc indicated by dots (from 
WC07) and squares (Leutenegger et al. 
2006). 



Waldron & Cassinelh (2007, WC07) analyzed fir 
line ratios for 17 OB stars, using plane-parallel Ku- 
rucz models, and found, in agreement with their ear- 
lier results, that high-Z ions are predominantly lo- 
cated closer to the photosphere. They highlighted 
a "near-star high-ion" problem. This problem can 
be illustrated by the Sxv hne in (^Puppis. WC07 
constrained i?fir for this line to < 1.22 i?*. This is 
in agreement with the results of Leutenegger et al. 
(2006) who estimated the minimum value of i?fir as 
l.lio 1 R*- The cmissivity of Sxv has its maximum 
at T wl6MK. The presence of plasma with such 
high temperatures at wl.li?*, i.e. close to the pho- 
tosphere, is contrary to the expectations from the 
shocked-wind model, but supports the base-corona 
model. 

WC07 notice a correlation between the radius 
where the ambient cool wind becomes transparent 
for the X-rays R{T^ind)~^, and the predominant lo- 
cation of He-like ion X-ray emission. This is illus- 
trated in Fig. 2. We use a spherically-symmetric 
PoWR atmosphere model of CPup to calculate 
-R(Twind)=l. The calculations are performed for two 
different M evaluated under the assumption of un- 
clumped wind (Repolust et al. 2004) and allowing 
for macroclumping (Oskinova et al. 2007). The mea- 
surements of by WC07 and Leutenegger et al. 
(2006) are also shown. As can be seen in Fig. 2, there 
is a good agreement between the measurements by 



WC07 of the /ir-inferred radii and i?(Twind)=l. How- 
ever, there is an apparent disagreement with the re- 
sults of Leutenegger et al. (2006). This can be due 
to a difference in measuring TZ and/or different ap- 
proaches to the analysis (see WC07). 

Because we preferentially see X-rays that are gen- 
erated close to the Twind = 1 surface (e.g. Ignace et 
al. 2000), the results of Leutenegger et al. imply a 
reduced mass-loss rate, while the WC07 results are 
consistent with the M estimates of Puis et al. (2006) 
and Oskinova et al. (2007). 

2.2 X-ray emission line profiles 
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Figure 3: O viii (thick blue) and Nex (thin red) hues 
in the spectrum of C Pup plotted versus 
wavelength in velocity units. The red- 
ward shoulder in Nc x line is due an un- 
resolved Fe emission feature. 

The analysis of fir-lines suggests that lines from 
lower Z ions arc formed further out in the wind. 
Therefore, assuming that the wind velocity mono- 
tonically increasing with the radius, lines from ions 
with lower Z should be broader than from ions with 
higher Z. Contradictory to this expectation, the 
observed widths of all X-ray lines are similar (c.g 
Kramer et al. 2003; Pollock 2007). To illustrate this, 
the lines of Ne x and O Viii in spectrum of ( Pup are 
shown in Fig. 3 (see also lines of C vi and Ne X in 
COri spectrum shown in Pollock 2007). Thus, there 
is a discrepancy between the location of the emission 
region inferred from the fir line ratios and from the 
fitting of the line profiles. 

Radiative transfer in homogeneous wind predicts 
that the redshiftcd emission from the back side of the 
wind is more suppressed by the continuum opacity 
along the line of sight, than the blucshifted emis- 
sion from the front hemisphere. This should lead 
to skewed, almost triangular-shaped profiles of the 
emission lines (MacFarlane etal. 1991; Ignace 2001; 
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Owocki & Cohen 2001). Because the atomic opacity 
increases at longer wavelengths, the lines of ions 
with lower Z are expected to be more skewed com- 
pare to the higher Z ions. However, the majority of 
observed stellar X-ray lines arc only slightly skewed 
and the shifts are similar for lines of all ions (WC07). 

There are a number of different hypotheses to ex- 
plain the puzzling properties of observed line pro- 
files. Pollock (2007) proposes that X-rays originate 
in coUisionless shocks. Two-component wind expan- 
sion is invoked by MuUan & Waldron (2006). Ig- 
nace & Galley (2002) show that if the X-ray emitting 
plasma is optically thick this would result in more 
symmetric line profiles (sec also Leutencgger ctal., 
these proceedings). Waldron & CassincUi (2001) 
suggest that nearly symmetric line profiles can be 
explained by low opacity of the stellar wind. 

Waldron & Cassinelli (2001) note that the wind 
opacity could be reduced if a large fraction of the 
wind were photoionized by X-rays. However, they 
observe that this is not compatible with the presence 
of low-ionisation species, e.g. C iv in the UV spectra. 
Our recent PoWR models of O star winds confirm 
this suggestion: there is no significant change in the 
ionisation structure of the wind when a realistic X- 
ray intensity is included. 

Waldron & CassincUi (2001) further suggest that 
if the adopted mass-loss rate is reduced this would 
lead to a better agreement between modeled and ob- 
served line shapes. This suggestion is supported by 
the modeling of line profiles emerging from a ho- 
mogeneous wind by Kramer et al. (2003) and Co- 
hen et al. (2006). They infer the wind optical depth 
by fitting a model with four free parameters to the 
lines observed in the X-ray spectra of C Pup and 
Ori. The optical depth in the smooth wind with 
a monotonic velocity law scales with a parameter 
n = K^M{vooR*)~^ (Owocki & Cohen 2001). Infer- 
ring r* from the model fits and assuming that the 
wind atomic opacity is a constant, Kramer et al. 
(2003) and Cohen et al. (2006) conclude that X-ray 
lines can be fitted only if M for CPup and COri 
are strongly reduced (see also Cohen, these proceed- 
ings). More realistic mass- absorption coefficients 
are presented in Oskinova et al. (2006). We use the 
best-fit values of r* obtained for different lines from 
Table 1 in Kramer et al. (2003) to estimate the mass- 
loss in CPup, however there is large disagreement 
when different lines are considered. 

The smooth-wind line profile formalism can be 
easily adapted for an inhomogeneous wind, where 
all clumps are optically thin at any given wavelength 
(so-called microclumping approximation). The em- 
pirical mass-loss rates based on diagnostics are 
lower when microclumping is adopted, compared to 
models that assume a smooth wind. The mass- 
loss rates obtained from the fitting of X-ray lines 
by Kramer et al. (2003) can be reconciled with the 
radio and Ha measurements only if clumping filling 



factors are very small (see Puis et al. 2006). Over- 
all, the use of microclumping approximation, at least 
in dense winds, is questionable (see Hamann, these 
proceedings). 

Waldron & Cassinelli (2001) also note that a 
clumped wind can be effectively optically thin even in 
stars with large M. This suggestion was exploited in 
Feldmeier et al. (2003) by waiving the microclump- 
ing approximation. 

Assume that the flow of clumps that constitutes 
the wind obeys the equation of continuity. The num- 
ber of clumps per unit volume is n{r)=nQv(r)~^r~'^ , 
where ng is a constant. The effective opacity of 
clumped wind, KcS = n{r)aciump'P, is the prod- 
uct of the clump number density n(r), the clump 
cross-section cTdump, and the probability of an X-ray 
photon being absorbed when it encounters a clump, 
T' = 1 - cxp (-r^^^^P) (Feldmeier et al. 2003). 

The optical depth across an average clump at 
the distance r is r^i™P = T^,R^Vaar~^n^^ . When 
^ciump ^ ^ -j-j^g effective opacity Koff does not de- 
pend on the atomic opacity k^, (because V — > oo). 
Therefore, the wind effective opacity becomes grey. 

To calculate emission line profiles, we assume for 
simplicity that the emission with emissivity rj orig- 
inates between some specific radii ri and r2 in the 
wind. The formal integral for the line profile reads 
Ffi ^ ry(r)r^e~'^"dr, where fi is the direction co- 
sine. 

Evaluating the wind optical depth as an integral 
over effective opacity along the line-of-sight z gives: 

Tw=no/ (1 -e~'"""""')d2. (1) 

To understand how the clump geometry alters the 
line profile shape, let us consider two extreme cases 
and compare isotropic clumps (balls) and anisotropic 
clumps in the form of infinitesimally thin shell- 
fragments (pancakes) oriented perpendicular to the 
radial direction, so-called a "Venetian blind" model. 
Assuming that the clumps keep constant solid an- 
gle as they propagate outwards, in the direction z a 
ball has the cross-section crocr^, while a pancake has 
(Toc|^|r^. Recalling that dz = dr/n and inserting 
the above expressions for the cross-sections in the 
Eq. (1), it is immediately clear that in the case of 
pancakes the integral Eq. (1) transforms into an in- 
tegral over r, while it stays an integral over z in case 
of the balls. The consequence is that pancakes yield 
nearly symmetric emission line profiles, similarly to 
the observed (Oskinova et al. 2006). 

Owocki & Cohen (2006) applied a porosity 
length formalism (Owocki et al. 2004) with = 
ft.(r)~^(l— exp(— T^'"™P)), where h{r) is a parameter. 
They studied a case of isotropic opacity, and con- 
cluded that large porosity lengths /i(r)wli?* are re- 
quired to reproduce the observed X-ray line profiles. 
However, Oskinova et al. (2006) computed synthetic 



4 



X-raying clumped stellar winds 



X-ray lines for anisotropic opacity. Stellar wind pa- 
rameters, i.e. K^, M, and v{r) were adopted from 
most recent stellar atmosphere models and it was as- 
sumed that the average separation between clumps 
at distance r is R<,v{r) jvoo, as compatible with pre- 
dictions of hydrodynamic simulations (e.g Feldmeier 
et al. 1997). Synthetic lines were computed without 
allowing any free parameter and compared with the 
observed lines. The remarkable similarity between 
synthetic and observed lines provides an evidence 
that the wind clumps are not optically thin, and are 
compressed in radial direction. 

3 X-ray emission of binaries as a 
probe of clumped stellar wind 

The collision of stellar winds in massive binary sys- 
tems (CWB) results in a bow shock that concaves 
around the star with the weaker wind. In close bina- 
ries the shocked material cools radiatively and hence 
the intrinsic X-ray luminosity is LxccMv'^ (Pittard 
& Stevens 2002). In the wide binaries with sufB- 
ciently large separation d, the shocked material cools 
adiabatically resulting in Lx'xAI^d~^v~^'^ (Stevens 
et al. 1992). 

Hydrodynamic studies of colliding clumped winds 
were conducted by Walder & Folini (2002) and Pit- 
tard (2007). Walder & Fohni considered WR140 
at periastron. They showed that the collision of 
clumped winds in a radiatively cooling system leads 
to the fragmentation of the colliding wind region 
(CWR). Carbon enriched dense WR-wind clumps 
with dimensions exceeding 10^^ cm can effectively 
cool and serve as the seeds for the formation of the 
dust. Pittard (2007) considered WR 140 at apas- 
tron when the CWR is likely to cool adiabatically. 
It was shown that in this case the CWR becomes 
highly turbulent. This leads to a rapid destruction 
of the clumps in the CWR. Therefore in adiabatic 
CWB, wind clumping does not affect the intrinsic 
X-ray luminosity. Therefore such systems may pro- 
vide a diagnostic of mass-loss that is independent on 
clumping. 

The X-rays generated in the bow shock propagate 
along the line of sight through the stellar wind of 
the foremost star, resulting in X-ray eclipses due to 
the orbital motion. The duration and the depth of 
the eclipses depends on the orbital parameters, the 
wind geometry, and its opacity. When the intrinsic 
X-ray luminosity from the bow shock is known, the 
stellar mass-loss rate can be probed by inferring the 
absorbing column from the X-ray spectrum. 

Schild et al. (2004) derived the mass-loss rate of 
the WC8 component of 7Vel in this way. The ob- 
tained M is a factor of four smaller than what is de- 
rived from spectral analysis with homogeneous wind 
models. This was interpreted in terms of WR wind 



clumping, with a volume filling factor of /v ~ 0.06. 
For ryCar, Pittard & Corcoran (2002) derived a 
mass-loss rate from an analysis of the X-ray spec- 
trum that is lower by a factor of few compared to the 
conventionally adopted value, thus indicating wind 
clumping. Similarly, Pollock et al. (2005) derived 
smaller than expected wind absorption from their 
analysis of the X-ray spectrum of WR140. 

Supergiant HMXBs consist of a supergiant OB 
star and (usually) a neutron star (NS), that orbits 
deep inside the stellar wind. An X-ray emission 
with a power-law spectrum results from the accre- 
tion of the stellar wind onto the NS. These X-rays 
photoionize the surrounding stellar wind. The wind 
X-ray spectrum shows a large variety of emission fea- 
tures, including a number of fluorescent lines. Sako 
et al. (2003) reviewed the spectroscopic results ob- 
tained with X-ray observatories for several wind-fed 
HMXBs. They concluded that the observed spec- 
tra can be explained only as originating in a wind 
where cool dense clumps are embedded in rarefied 
photoionizcd gas. Sako et al. (1999) constrained the 
volume filling factor of clumps in the wind of Vela 
X-1 as /v «0.04. Similar conclusions were reached 
by Van der Meer et al. (2005). They studied the 
X-ray light curve of 4U 1700-37 and found that the 
feeding of the NS by stellar wind clumps explains 
the observed stochastic variability. The clump sepa- 
ration at the distance 2R^, was constrained to 0.4i?*. 
Recently, the advances in the 7-ray and X-ray obser- 
vations lead to the discovery of new HMXBs that are 
highly absorbed and display fast stochastic variabil- 
ity. This spectral and temporal variability is con- 
sistent with strong wind clumping. The inferred 
clumping parameters are similar to those obtained 
from the analysis of X-ray emission in single stars 
and in CWBs (see Romero et al. and Walter et al., 
these proceedings). 



4 Conclusions 



The spectral and temporal properties of the X-ray 
emission from single stars, colliding wind binaries 
and high-mass X-ray binaries provide consistent ev- 
idence of wind clumping. The multitude of data 
requires strongly clumped winds where clumps are 
separated by a few tenths of the stellar radius in 
the wind acceleration zone. The approximation that 
clumps are optically thin at all wavelengths in the X- 
ray band is not justified and cannot be applied uni- 
versally. Clumped stellar winds have a reduced opac- 
ity for X-rays, compared to a homogeneous wind, 
due to a reduction of empirically determined mass- 
loss rates and to the effect of wind porosity. 
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